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Abstract 

The effects of classical gravitational radiation in models with large extra dimensions 
are investigated for ultra high energy cosmic rays (CRs). The cross sections are 
implemented into a simulation package (SENECA) for high energy hadron induced 
CR air showers. We predict that gravitational radiation from quasi-elastic scattering 
could be observed at incident CR energies above 10^ GeV for a setting with more 
than two extra dimensions. It is further shown that this gravitational energy loss 
can alter the energy reconstruction for CR energies -EcR ^ 5 • 10^ GeV. 
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1 Motivation 



One of the major problems in modern physics is to combine quantum physics 
and gravitation. The most promising candidate to solve this problem seems 
to be string theory because it includes all symmetry groups of the stan- 
dard model, however, it introduces additional space dimensions [1]. From the 
present (non-)observation of these dimensions it is concluded that they are 
compactified and their size was assumed to be on the Planck scale. However, 
string theory itself does not give any stringent constraints on the size of the 
extra dimensions other than non-observability so far. 
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This opens a possibility for gravity-only extra dimensions, assuming that only 
gravitons arc able to enter these additional dimensions. The size of these di- 
mensions is then primarily constraint by direct measurements of the gravita- 
tional inverse square law. However, the strength of gravitational interactions 
has only been measured down to a scale of some micrometers [2,3]. Below 
this scale some modification of the gravitational force law is still possible. 
This allows to introduce a new fundamental scale of gravity Mf <^ Mpianck 
in the TeV-range. These Large Extra Dimensions models (LXDs) have com- 
pactification radii up to ~ /xm and a new fundamental scale in the TeV range 
[4,5,6]. If the new scale of gravity is indeed not far beyond the electro-weak 
scale, it might be reachable in future collider experiments at LHC and be- 
yond or in cosmic ray observatories like AUGER or ICECUBE. At center of 
mass energies above the new fundamental scale Mf the LXDs become im- 
portant and first observable effects of quantum gravity might be observable. 
A multitude of new effects have been predicted in the recent years to ob- 
tain deeper insight into this exciting field ranging from Kaluza-Klein graviton 
production and increased neutrino cross sections to black hole production 
[7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27]. 

Apart from these promising phenomena, present day observations also allow 
for direct and indirect constraints on the new fundamental scale and the size 
of the LXDs: 

• A fundamental scale in the TeV-range does not allow a single LXD as the 
compactification radius would get in the order of the size of our solar system 
and we would have noticed. 

• For higher numbers of LXDs the compactification radius does not exceed 
the /im-scale as shown by the mentioned short distance experiments. 

• Even smaller distances can be probed in collider experiments. At the mo- 
ment the highest center of mass energies are reached in hadron-hadron col- 
lisions at Tevatron with 2 TeV cm-energy which sets a bound of Mf > 1.0 
TeV for 2 down to 0.7 TeV for 6 LXDs [28,29]. 

• Other bounds are set by supernova explosions [30,31,32,33]. The cooling of 
the supernovae is modified by the production of gravitons which is very 
sensitive to the number of LXDs in the reached energy region. For 2 LXDs 
a bound of about 500 TeV is obtained. However, for any higher number of 
LXDs the bounds are lower than from the collider experiments. 

In the last years it was recognized that cosmic rays (CRs) provide an excel- 
lent laboratory to study the onset of physics beyond the standard model. The 
initial energies of UHECRs exceeds 10^^ GeV ~ 100 - 1000 TeV) 

and might therefore allow one to probe TeV scale large extra dimensions 
[34,35,36,37,38,39,40,41,42,43,44,45]. Especially the radiation of gravitational 
waves in models with large extra dimensions is strongly enhanced compared 
to standard general relativity. As will be discussed later, this might lead to 
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observable signatures and modifications of the reconstructed flux and energy 
estimates for UHECRs. 

To explore the effect of gravitational energy loss in cosmic ray air-showers, we 
use the differential cross section for gravitationally radiated energy in quasi- 
elastic scatterings as calculated in [46,47] for 2 ^ 2 processes. Then we apply 
the result to cosmic ray air shower simulations and extract the impact on 
cosmic ray observables. 



2 Gravitational radiation from quasi-elastic scattering with extra 
dimensions 



First estimates to study effects of gravitational energy loss of CRs due to the 
presence of extra dimensions were explored by [38] . There, the presence of large 
extra dimensions was incorporated into the well known results from general 
relativity [48] by a change of the phase space seen by the emitted gravitational 
wave. The additional phase space factor for the emitted gravitational wave was 
given by 

ddK'^d) - dV{d/2)T:<iin<i-^- ^ ^ 

Note that Qq — I. Where R is the compactification radius of the extra dimen- 
sions in the ADD scenario given by [4] : 

d+2 2 

R^Mf MpV (2) 

Here, Mf is the new fundamental scale and Mpi is the four dimensional Planck 
mass related to the gravitational constant by Gn = V-^pi- 

This method led to a strong modification of the reconstructed energy spectrum 
and the authors concluded that the steepening of the CR energy spectrum 
around 10^^'^ eV (the "knee") might be due to gravitational energy loss. 

However, from our present studies (see also [46]) it seems that a calculation of 
the effects of the gravitational energy loss requires a more elaborate treatment 
as will be discussed now. The simplified treatment can be improved by direct 
calculation ^ of the gravitational energy loss in a A'" — > M scattering process 



^ We will use the following notations: 4 + d space-time vectors will he x = {xq^x)^ 
where the spatial part can be split again into a three dimensional and a d dimensional 
part X = (cc, x_[_). 
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as given by [46,47]: 



(3) 

The P(7v) are the momenta of the coUiding particles and the factors rj^ are 
defined by 

— 1 for initial state particles. 
W = \ (4) 
+1 for final state particles. 

Thus, in the case of a 2 — > 2 collision the index runs from 1 to 4. Before 
we continue, we want to point out that equation (3) follows from classical 
considerations and is not derived from any form of quantum theory of gravity 
(e.g. loop quantum gravity, SUGRA or string theory). However, we believe 
that it can account - at least semi-quantitatively - for the major effects of the 
gravitational energy loss. 

Next, we integrate Eq. (3) with the help of the Mathcmatica package FcynCalc 
[49,50]. Difficulties for the dl^s+d integration arise from the P ■ k terms in the 
numerator. The protons are bound to the brane and the product P ■ k gives 
for example for one of the incoming protons: 



Pi • A; = P^ko - Pik - = |p|A;o 



1 + 



\ \P'\ \ kl 



k'^ 

rfcos0fe . (5) 



For kl/k^ and (pk ^ 0, Pi ■ k becomes small and the denominator in Eq. 
(3) is only regularized by mp/|p^|. 

We introduce the Mandelstam variables s, t and u by 

S^{Pl+P2f, t^iPi-Psf, U^{Pi-P,f. (6) 

It is convenient to perform a coordinate transformation to rewrite Eq. (3) in 
terms of spherical coordinates in three dimensional space and the d extra- 
dimensional coordinates separately: 



dE k'^+^dE k'^+HE {k" + klY'^+^^l^dE 



dkodVL^+d dk^+<^ dk<^dk^ k^dkdQ.z\kd\'^-^dkddVLc. 
Solving this for the new integration variables yields 



(7) 



dE dE k^kt^ 



dkdQsdkddQd dkodQs+d {k^ + kj) (<^+2)/2 ' 



(8) 



The first term on the right side can be approximated by Eq. (3) as soon as 
the wavelength of the gravitational wave is smaller than the compactification 
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radius R of the extra dimensions and the gravitational wave can propagate 
freely into the bulk. Rephrased as a condition for this constrained becomes 



1^1 > Mf 



Ml. 



(9) 



A lower bound on |A';| is not relevant for the energy loss discussion, because the 
major contribution to the radiated energy comes from the high energy (i.e. 
large \k\) part. To calculate the energy loss due to the gravitational wave emis- 
sion one has to perform the dfl^ = sin (pkzdfpkdfpkzi the dVLdi the dk and the dk^ 
integrals. However, the rather steep t dependence of the elastic Proton-Proton 
cross section allows us to simplify these integrals, because the physically rele- 
vant processes are dominated by small |t| < contributions, with irip being 
the mass of the Proton. Thus, one can expand Eq. (3) for small \t\. This gives 
for the part E/ j ■ ■ ■ in Eq. (3) containing the sums over external momenta 



i,j 



= -St{{kl - kl)\Aml - s)s'cos{ct>kf 

+ kl{kl - kl)s{Aml + s)cos(0fc)^ [A;^s(-8mJ - Amis + s^) 

- 4A;g(32m^ - 14mjs - 3m^s2 + 

+ {kl - kj)s{-8m^ - Amis + s^) cos(20fc) 

- 0.5 [{kl - klfs^cos{(Pkf {klsiSml - Am^ + s^) 
+ kl{-128ml + 88mp + I2mls^ - 7s^) 

+ {kl - kl)s{8ml - Amis + s^) cos(20fe))" 



- 0.5 



kt{Aml + sf {kls{8ml - Amis + s^) 



+ kl{-l28ml + 2Amis + Umls"^ - 3s' 



+ 



{kl - kl)s{8ml - Amis + s^) cos(20fc))] } 



/ \kl{-Amp^^s) 



Aml + s+[-l + f^ 



S COs{(f)k 



n4 



(10) 



For ~ the radiation does not propagate into the extra dimensions and Eq. 
(10) reduces to the well known classical limit. From Eq. (10) one can see that 
for k'^/kls > Ami regularising part in the denominator is not rrip/s any 
more and a Taylor expansion of Eq. (10) around m^/s = is allowed. This 
expansion has a large validity region for ultra high energy collisions because 
it just demands that 



2 s 



fill 



This approximation also fulfils the condition in Eq. (9). After performing the 
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integration ^0^2, this series gives 



dE _ t 2ki-^k^[{kl- kl) cos {2^k)~kl 



dkddkd(t)kd^d {2uY+^Mj+^ [kl + {kl - kl) cos {<Pk)Y 
Next we perform the 4>k integration, 

dE _ t A;^-'fe'(2A;2 + 3fc') 



(12) 



(13) 



dkddkdQd {2TTYMf+'' (A;i + fe')2 
The integration over the d dimensional unit sphere fid gives a factor 27r'^/^/r((i/2). 

dE ^ t kt'e{2kl + 3e) 

dkddk 2d-^7r'^/^r{d/2)Mf+^ (kd + ky ' ^ ^ 

Next, the kd and the |fc| integration can be performed with respect to the 
integration hmits k^ + k'^ < kmax and Eq. (11). This calculation can be done 
explicitly for two, four and six extra dimensions: 

E{t, d = 2) = -kl,J [5V2 - log (1 + v^)] /iUirMf), 

E{t, d = 4) = -kl,J [V2 - log (1 + v^)] /(levr^M^), (15) 

E{t, d = 6) = -kl^^t [11V2 - ISlog (1 + V2)] /(17927r3M«). 

Let us now discuss the relation between this result and those obtained in 
earlier publications: 

• In Ref. [46] the gravitational wave was assumed to have a momentum vector 
only in the direction out of the brane, thus the denominator in Eq. (3) sim- 
phfies to Pjk — Ejko- After integrating over ko (which is not strictly correct, 
because the problem is not spherically symmetric in 3+d spatial dimensions 
any more) the result shows the same t and kmax dependence as Eq. (15). 
The different factors are due to the simplification in the integration. 

• The phase space argument used in Ref. [38] leads to the same kmax de- 
pendence. However, the pre-f actors differ and even more striking the result 
derived in [38] has no t dependence (but s instead). Therefore, this ap- 
proach leads to drastic overestimation of the gravitational energy loss in 
high energy cosmic rays, as we will see in the following sections. 



3 QuEisi-Elastic hadr on- nucleus scattering 



In order to calculate the energy loss due to gravitational wave emission in air 
showers at high energies one has to know the elastic scattering cross section 
d(Teiastic/di. We construct it from the hadron-nucleon scattering cross section. 
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In the impact parameter representation, the profile function is defined as 



T{s,b) 



l^total 



exp 



2B. 



elastic , 



(16) 



2 27r-Belastic 

where B is the elastic scattering slope and o"totai is the total hadron-nucleon 
cross section. The profile function is related to the amplitude, 



r{s,b) 



Stt 



A{s,b) . 



(17) 



The profile function Eq. (16) can be expressed via the eikonal [51] function x 

by 

r(s,6) = l-exp[zx(6)] , (18) 
and is related to the phase shift of the scattered wave. 

The Fourier transform relates impact parameter space to momentum space by 



A{s, t) = ^ J d^bA{s, b) exp(-i6qf) , 



(19) 



with t = —q^ ~ —Q^- The amplitude gives the differential elastic cross section: 

dCelastic 



dt 



167rs 



dt 



exp(5eiastic^) ■ 



(20) 



t=Q 



Prom the last expression we can define the elastic scattering slope 



B = 



d A dfJclastic 

di \ 



■A 

dt J 



(21) 



t=o 



Inserting Eqs. (19) and (20) into Eq. (21), the scattering slope becomes 
B 



2^1 / d^bA{s,b)e-''"i\ 

I / d2bA(s,b)c-''"3| 



t=0 



-2qdq 



I / d26yl(5,b)(l+(-i6q)+i(-i6q)2-...)| 



I / d^bA{s,b)\ 
I J d'^bA(s,b){b'^ cos'^ <p)\ _ J dbb^r{s,b) 
I / d^bA{s,b}\ ~ 2 J dbbr(s,b) 



t=0 



(22) 



where we have expanded the exponential and kept only the third term. The 
first term does not depend on q and the second term in the expansion vanishes 
due to symmetry. 

To obtain the scattering slope of a hadron-nucleus collision we replace the 
hadron-nucleon scattering profile function by 



^hA{s, 6) = 1 - exp 



-tJ2xi{s,bi) 



i=l 



1 - 



1-Tk(6) 



(23) 
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10^ 10^ 10^^ 10^ 10® 10^ 

Invariant Energy [GeV] 



Fig. 1. The hadron- nucleus slopes (thin lines) and the hadron-nucleon slopes (thick 
lines) as a function of the collision energy in the center of mass frame. Primary 
particles are Protons (full lines), vrs (dashed lines), and Kaons (dotted lines). 

where is obtained from the Glauber-Gribov formalism [52,53,54] by con- 
volution of the thickness function with the hadron-nucleon profile, 

f^{b) = jd'cT^{bnb-c) , ^^^^ 

TA{b) = \Jdzpiz,b). 

Using Eq. (24), the scattering slope for elastic hadron- nucleus collisions Eq. 
(22) is calculated as a function of the collision energy. The underlying hadron- 
nucleon scattering slopes are taken from the SIBYLL model [55,56]. Fig. (1) 
depicts the underlying hadron-nucleon slopes (thin lines) and the calculated 
hadron-nucleus slopes (thick lines). The hadron- nucleus slopes are clearly 
higher than the hadron-nucleon slopes at the same energy. However, the ratios 
of the two slopes decreases with increasing energy. 



4 Gravitational radiation from high energy cosmic rays 

After the derivation of the basic equations in the previous sections, we are 
now ready to calculate the amount of energy that is emitted into gravitational 
radiation by a high energy proton propagating through the atmosphere. 
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unitarity-bound - - - 

6=2 

d=4 

d=6 



"10^ 10^ 10® 10^° 10" 

Energy [GeV] 

Fig. 2. Energy loss (in GeV/(g/cm^)) of a proton propagating through the atmo- 
sphere as a function of the lab-frame energy for Mf = 1 TeV and d = 2,4, 6. 

The differential energy loss is given by 



10' 



10"' 



dx 



s, d) 



/ dt 




dt 



E{t,s,d) 



da^ 

\ J dt^ 



(25) 



where A is the mean free path for elastic scattering of the projectile in units of 
g/cm^ and da^^/dt is the differential hadron-nucleus cross section. For cosmic 
ray calculations it is convenient to calculate the energy loss E in the laboratory 
frame. The corresponding Lorentz transformations are given in App. A. 

Figs. 2 (for Mf = 1 TeV) and 3 (for Mf = 2 TeV) show the differential energy 
loss of a Proton propagating through the atmosphere as a function of the 
initial energy in the laboratory frame. The short dashed, dotted and full lines 
give the results for two, four and six extra dimensions, the long dashed lines 
show the unitarity bounds. For large initial energies, a higher number of extra 
dimensions leads to an enhancement of the gravitational energy loss. However, 
with increasing fundamental scale Mf the effect is much weaker as shown in 
Fig. 3. Note, that the result is cut-off dependent as fcmax is not determined 
from first principles. For the present study, we have chosen femax = 
which is the maximal value consistent with energy conservation in the picture 
of a gravitational wave being emitted by one of the outgoing states. The 
comparison of these results with [38] shows that an approximation of the 
effects of extra dimension with a simple phase space argument does yield a 
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Energy [GeV] 



Fig. 3. Energy loss (in GeV/(g/cm^)) of a proton propagating through the atmo- 
sphere as a function of the lab-frame energy for Mf = 2 TeV and d = 2,4, 6. 

similar shape for the energy loss as those shown in Figs. 2 and 3. However, 
the omission of the correct kinematics of the energy loss Eq. (15) results in a 
dramatic overestimation of the gravitational energy loss effect by several orders 
of magnitude. In addition, the simple extension of the standard formula with 
a modified phase space factor on the integrated cross sections results in a 
violation of the unitarity bound. 

Even though the energy loss into gravitational waves in our (very optimistic) 
scenario is much lower than expected from previous approximations, it might 
still have observable consequences for very high energy cosmic rays. Therefore 
we implemented Eq. (15) and the elastic cross sections into a complete cosmic 
ray air shower simulation (SENECA) [57,58] to study the modifications of the 
shower properties in detail. 

Fig. 4 gives the relative energy loss as a function of the incident energy E. 
The calculation is averaged over incident zenith angles dcos(6') in the range 
0° < ^ < 60°. The full lines indicate the calculations with six extra dimensions, 
while the dotted lines show the results for four extra dimensions [Mf = 1 TeV 
is shown by thick lines, Mf = 2 TeV is shown by thin lines). For the case of two 
extra dimensions, deviations from the non-modified shower properties are very 
small even for the most optimistic cases. However, for four extra dimensions 
first deviations from the standard calculation become visible at energies higher 
than 5 ■ 10^° GeV. For d = 6 the gravitational radiation becomes sizeable and 
already leads to deviations around 5 ■ 10^ GeV. At the highest energies, the 
integrated relative energy loss due to gravitational radiation might even exceed 
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0.25 



,12 



Energy [GeV] 



Fig. 4. Relative energy loss into gravitational radiation as a function of the incident 
cosmic ray energy E for d = (4,6) and Mf = (1, 2) TeV. 

20% of the initial particle energy. 

In present day experiments, e.g. AUGER, this gravitational energy loss would 
show up as a decrease in the number of observed secondary particles. The 
multiplicity of secondary particles Nsec{E,x) is directly observable in fluores- 
cence experiments and is a key observable to estimate the cosmic ray's initial 
energy. Any non-visible energy emission results in an underestimation of the 
initial energy in the energy reconstruction procedure. Thus, it has an impact 
on the interpretation of the measured cosmic ray flux in dependence of the 
incoming particle energy. 

How big is the distortion of the reconstructed flux due to graviton emission 
quantitatively? Neglecting fluctuations, for a given incoming flux F = dN/dE, 
the measured flux F' = dN'/dE' depends on the reconstructed energy E'{E). 
By identifying the integrated fluxes N'{E') = N{E) one finds 



For an incoming flux F = kE~^ the flux reconstruction is shown in Fig. 5. In 
all scenarios {d > 4) gravitational wave emission might indeed influence the 
energy reconstruction above 5 ■ 10^ GeV. For ultra high energy cosmic rays 
even an apparent cut-off seems possible ^ , because the relative amount of non- 

^ Note the linear scale on the y-axis, thus the suppression does not mimic the GZK 



F'{E') 



dN' dN (E) dE 



F{E). 



dE 
dE' 



(26) 



dE' dE ' dE' 
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1.05 




Energy [GeV] 

Fig. 5. Reconstructed flux F'E'^ /k as a function of the incident cosmic ray energy 
E' for d = (4, 6) and Mf = (1, 2) TeV. 

visible energy increases strongly with increasing energy. Hence, for UHECRs 
the interpretation of experimental data might have to be modified in scenarios 
with large extra dimensions. 

Presently available data from Hires and AGASA do not allow one to observe 
the predicted suppression pattern, because even in our most optimistic sce- 
nario the flux is reduced only by a factor of 0.5 for the highest energies. How- 
ever, with the expected high statistics data from the Pierre Auger Observatory 
a detailed exploration of this phenomenon might be possible. 

As a remark, we want to point out that in our calculation, gravitational wave 
emission does not give new insights into phenomena at lower energies {E < 
10^^ eV) and can not be considered as a candidate to explain the famous knee 
in the cosmic ray spectrum. 



5 Conclusion 

The energy loss into gravitational waves is calculated for ultra high energy 
cosmic rays. In contrast to previous estimates, quasi-elastic particle scattering 
in the ADD scenario with 4 or 6 extra dimensions has no observable influence 
on the properties of cosmic ray air showers at incident energies below 5-10^^ eV. 

cut-off. 
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Thus, the emission of gravitational radiation can not be used to explain the 
steepening of the cosmic ray spectrum at the "knee" {E ~ lO^''''^ eV). For two 
large extra dimensions, the studied effects are generally too small to lead to 
any observable effect. 

However, for energies above 5 • 10^^ eV and Mj < 2 TeV, d>A gravitational 
energy loss during the air shower evolution can be sizeable. This might result 
in an underestimation of the reconstructed energy for ultra high energy cosmic 
rays as studied by Hires, AGASA and the Pierre Auger Observatory. 
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A Energy loss in the lab system 



Equation (15) provides the gravitationally radiated energy in the centre of 
mass frame of the reaction. To transform the kinematic variables to the lab- 
oratory frame with a target Proton at rest one has to apply the Lorentz 
transformation matrix 



A 



^ Cosh{ri) Smh{ri) \ f ^ yTTe^ 
^ Sinh(7y) Cosh(7y) 




2mp 

VL 

2mr, 



(A.l) 



which acts on the t and the z (i.e. longitudinal) component of the A + d 
dimensional vector. All the other (transverse) components remain unchanged. 
Eq. (15) gives the energy E and momentum k of the gravitational radiation 
emitted from one of the interacting particles (pi, p2). For different momentum 
directions k/\k\ the Lorentz transformation Eq. (A.l) gives different energy 
losses in the lab-frame. To avoid this complication we use a mean value of the 
left over four momentum p' of the scattering particles. If the energy is radiated 
away from particle i we define p'i — Pi — k. Averaging over these cases yields 
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Using the symmetry of a 2 — > 2 scattering in the centre of mass system we 
find 

FcM= (v^-^0,0,0,0,...). (A.3) 

Because p'cu has no z component, the mean left over energy in the laboratory 
system becomes 

p'iab = A-FcM= (^(v^-M,0,0...^. (A.4) 
From Eq. (A.4) the mean energy loss in the lab system is obtained as 
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